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ABSTRACT 

The influence of an external magnetic field on a quasi one-dimensional system with a 
charge density wave (CDW) instability is treated within the random phase approximation 
(RPA) which includes both CDW and spin density wave (SDW) correlations. We show that 
the CDW is sensitive to both orbital and Pauli effects of the field. In the case of perfect 
nesting, the critical temperature decreases monotonously with the field, and the wave vector 
of the instability starts to shift above some critical value of magnetic field. Depending 
on the ratio between the spin and charge coupling constants and on the direction of the 
applied magnetic field, the wave vector shift is either parallel {CDWx order) or perpendicular 
{CDWy order) to the most conducting direction. The CDWx order is a field dependent linear 
combination of the charge and spin density waves and is sensible only to the Pauli effect. The 
wave vector shift in CDWy depends on the interchain coupling, but the critical temperature 
does not. This order is affected by the confinement of the electronic orbits. By increasing 
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the relative strength of the orbital effect with respect to the Pauli effect, one can destroy 
the CDWy, establishing either a CDWx, or a CDWq (corresponding to perfect nesting wave 
vector). We also show that by increasing the imperfect nesting parameter, one passes from 
the regime where the critical temperature decreases with the field to the regime where it is 
initially enhanced by the orbital effect and eventually suppressed by the Pauli effect. For a 
bad nesting, the quantized phases of the field-induced CDW appear. 
PACS numbers: 71.45L, 71.70E, 75.30F 
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1 Introduction 

The open and almost flat Fermi surface that characterizes the quasi-one-dimensional 
(QID) electronic systems gives rise to the formation of charge (or spin) density waves 
|1|, 0, ^. Moreover, the external magnetic field couples to the spin (via Pauli term) 
and to the orbits (via Peierls substitution in the Hamiltonian) of the electrons. This 
coupling affects the properties related to density wave [DW) formation like the order 
parameter, the critical temperature and the wave vector of instability. The scale for 
the Pauli impact in the momentum space is the wave number qp = ^bH/vf, while the 
orbital effect enters through the inverse magnetic length Qo = ebH cos 6, where 6 is the 
inclination of the magnetic field H from the transverse c-direction in the {b, c) plane (a 
plane perpendicular to the chains). The ratio of these two characteristic wave numbers 
T] = qo/qp = ebvp cosO / hb is of the order of unity in real materials. It will play an 
important role in the phase diagram for the CDW in a magnetic field. 

The Pauli term introduces a finite coupling between the CDW and the component 
of the SDW parallel to H, and may lead to a finite, magnetic field dependent, shift 
in the wave vector of instability It is therefore necessary to treat CDW and 
SDW together. A simple relevant model is the extended Hubbard or ((71,(72) model 
with coupling constants Us = (72/2 and Uc = {2gi — (72)/2 for the SDW and 
CDW respectively. Since the Pauli term mixes the CDW with the SDW, the ratio 
u = —Us/Uc will be the second relevant parameter for the CDW phase diagram. 

The Pauli term breaks the rotational symmetry of the complex vectorial SDW 
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order parameter, constraining its direction perpendicularly to magnetic field. With this 
constraint taken into account, the SDW phase diagram depends only on the orbital 
coupling, provided that the system is perfectly magnetically isotropic in the absence 
of magnetic field. However, the fluctuations of the component of SDW parallel to H 
around its zero value remain affected by both, Pauli and orbital coupling. Moreover the 
Pauli term introduces a finite coupling between these fluctuations and the noncritical 
CDW fluctuations. 

The influence of a magnetic field on the CDW systems is even richer, because 
both Pauli and orbital effects can affect the CDW ordering. This fact is of direct 
experimental interest, since e. g. the critical temperature can easily be measured. 
Furthermore, there is a finite magnetic field at which the wave vector of ordering starts 
to vary with the magnetic field. The description of these features, together with the 
interesting CDW — SDW mixing, is the main objective of the present detailed analysis. 

The various aspects of the interaction between the electrons in QID systems and the 
external magnetic field have been already subjects of numerous analyses. The quadratic 
decrease of the mean-field critical temperature in one- dimensional CDW systems due 
to the Zeeman splitting was proposed theoretically 0, and found experimentally in 
the organic compound TTF — TCNQ The recent very precise measurements 
in Per2[Au{'mnt)2] show the decrease of which differs considerably from the 
theoretical value 0. The effect of the Pauli coupling on the CDW order parameter can 
be formulated as a breaking of degeneracy of two density waves, those with parallel and 
antiparallel spin with respect to H, each component being a CDW-SDW hybrid. This 
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is reminiscent of the treatment of two coexisting CDWs with overlapping electronic 
bands P, ^, |l^. The coupling of two CDWs with different wave vectors may stabilize 



a soliton lattice in the relative phase of two waves |T0 |. 



On the other side, the orbital coupling alone leads to an increase of the critical 



temperature for CDWs pA\ , [T^, |T^]. Such an increase was observed in e. g. NbSes [|14 . 
The aim of the present work is to introduce both, orbital and Pauli couplings, into the 
RPA calculation of the DW matrix susceptibility, and to determine some mean-field 
properties, in particular the phase diagram for CDW systems in a magnetic field. 

In Sec. 2 we derive the RPA results for DW response functions in the form of a 
4x4 matrix. In Sec. 3 we analyse in detail the phase diagram for CDWs in the case 
of a perfectly nested Fermi surface. In particular we consider the influence of the 
parameters u, t] and of the interchain hopping on the critical temperature, the wave 
vector of the instability, and the CDW-SDW coupling. We also shortly discuss the 
effects of the imperfect nesting on the critical temperature as a function of magnetic 
field. The concluding remarks are given in Sec. 4. 
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2 Model 



Quasi-one-dimensional electrons in an external magnetic field are usually modelled 
by the anisotropic two-dimensional Hamiltonian 

h 



Ho = — y dqy J dx *^(a;, qy) [Hm + HQiD,orb + Hpaun] "^{x, Qy) 



with 



Hid = ivpPsdx (2a) 
HQiD,orb = '^hPs sm{qyb - Qox) + 2^ cos 2{qyb - qgx) (26) 
Hpauii = —asiisH (2c) 

Here and are four-component fermion fields, 

where the indices t, J, span the spin space and (7j are corresponding Pauli matrices. 
Indices +(— ) denote the right (left) Fermi surface with the states defined with respect 
to ±Q/2, where Q = {2kF, Tr/b) is the wave vector of perfect nesting reahzed for = 0, 
and Pi's are the Pauli matrices in that space. The chains lie in the xy plane and are 
parallel to the x axis, b is the latice constant in the y direction. The longitudinal elec- 
tronic dispersion given by Hid is linearized in the vicinity of the Fermi wave numbers 
±kF , with Vf being the longitudinal Fermi velocity, tb is the hopping integral between 
nearest neighboring chains and parametrizes the imperfect nesting. The spin space 
is chosen to have the third component parallel to H. 
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Let us now introduce the relevant interaction part of the Hamiltonian. Since 

the further considerations are hmited to the 2kp RPA response, it is sufficient to 

keep only the contributions with bilinearly coupled electron-hole operators for spin and 

charge density waves. They are given by 

^int = fdxJ2 [-UsM\R) ■ M(R) + [/eMl(R)M4(R)]. (3) 

The two-fermion operators in eq.(|^) are defined by 

Mi = ^V+^i^ ^ = 1,2,3,4 (4) 

where (T4 = /. The first three components {i = 1, 2, 3) define the complex SDW 
vector amplitude M, while the fourth component M4 is the complex CDW scalar 
amplitude. The SDW and CDW coupling constants in eq.(|^) are related to the usual 
backward {gi) and forward ((72) electron-electron coupling constants by f/^ = (72/2 and 
Uc = {2gi — (72)72. We shall specify later the range of these constants for the most 
interesting physical cases relevant for our analysis. 

The mean field (MF) critical temperature for the spin or charge density wave is 
defined as the temperature at which the corresponding RPA susceptibility diverges. In 
our case the Pauli term introduces a finite coupling between the component of SDW 
parallel to the magnetic field (M3) and the CDW (M4). This coupling is appropriately 
treated by introducing the DW susceptibility matrix, with the elements defined as 
retarded correlators 

X.,(q,t-t') = (M.Mj) = -^^(t-0([M,(q,t),Mj(q,0]) , ^,J = l,-,4 , (5) 
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where q is the deviation of the wave vector from (2fcir,7r/6). The RPA result for this 



matrix is 



[xii(q>^)] 



g Xo(q.'^) 













with 














/II 



x,{VTTs^-UsXg) 



(6) 



Xg = V^T(q,w)xi(q,^) 



XTa(q>^n) = Xoiqx±2qp,qy, 



s = [xt(q>^) - Xi(q,^)]/2xg 



/II 



= l + iU,-Us)XgVl + 6^-U,Usxl 
= l-UsXo{<i,uj) . 



(7a) 
(76) 
(7c) 

(7e) 



Xo(q, 1^) is the susceptibihty which includes orbital contributions of a magnetic field, 

oo 

Xo(q,cu„) = J2 Piqx-lqo,(^n)ifiqy), (8) 

l=—oo 

where P{k — Iqo, Wn) is the one-dimensional bubble. The coefficients Ii{qy) bring in the 
orbital quantization due to the finite transverse dispersion (2b) [jl2|, |15|, 



^i{qy) = 1^ Ji-^i' sm -^\ Jv \ cosgj^ft 

^ \VFqo 2 ) \VFqo , 



(9) 



where J; are Bessel functions. 



Even without further diagonalization of the matrix (|^), it is evident that the critical 



temperatures for the condensation of density waves follow from the conditions 

r(qx,T,^) = (10) 

and 

/"(qii,T|) = o, (11) 

where all functions have to be taken in the static {u = 0) limit. is the critical 
temperature for the SDW with the orientation of the spin perpendicular to H, i.e. for 
the degenerate components Mi and M2. TJ is the critical temperature for the hybrid 
of the CDW and the SDW with the spin parallel to H, i.e. of the coupled block 
(M3,M4) in the matrix [xjj(q, = 0)]. The corresponding wave vectors q_L and qy of 
the ordering are those which maximize the respective critical temperatures and Tj' . 
The true critical temperature of the DW instability is equal to max{T^ ^ TJ'}. 

Having in mind real systems, it is appropriate to distinguish the most important 
situations realized for two characteristic interaction schemes. In the case of repulsive 
interactions {Us > , f/c > 0), usually analysed in terms of the Hubbard model 
{Us = [/c > 0), the stable ordering following from (^ is the SDW one, determined 
by the condition (p^O]). In other words, as far as the system possesses the internal 
magnetic isotropy, there is no effect of Pauli coupling on the ordering. Its spin is 
oriented perpendicularly to H, while the wave vector is given by q_L = in the case 
of the good nesting {t'^^ <^ Tc), and may pass through the well-known cascade of phase 
transitions due to the orbital effects when the deviation from the good nesting is large 
enough (f, > Te) Pi H, m. 



In the case of predominant electron-phonon interaction {Uc < 0, \ Uc\ > Us > 0) 
the system prefers the CDW ordering. As it is obvious from eq.(|]), the off-diagonal 
matrix elements vanish in the absence of a magnetic field. The instabihty condition 



(|TT|) then reduces to 1 + f/cXo = 0, and the ordering involves only the CDW component 
M4. For finite magnetic fields the relation (|ll|) contains contributions originating from 
both, orbital and Pauli terms in the Hamiltonian (|1],0). The former enters through the 
bubble susceptibility (||), while the latter introduces the CDW - SDW (i.e. M4 - M3) 
hybridization measured through the parameter 5. As it is seen from eq.(0), 5 is finite 
if g^; 7^ 0. More explicitly, after diagonalizing the M3 — M4 block of the matrix (^ the 
normal components of the "vector" (^ read 

M_ = 1[5M3 + AM4] 

M+ = 1 - AM3] , (12) 

with 



N = ^5^ + A2 

and 

^"-2 vv 2 

while the corresponding diagonal susceptibilities are 



(13) 



In these equations we have defined 



U = -U, , v = UJU = -UJU, (14) 
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as a convenient parametrization of coupling constants for the problem of the 

CDW in the magnetic field. The value of v depends on the interactions which 
participate in the Hamiltonian (|l|, ^). The global phase diagram [jl| at = in 
(z/, f/)— space is shown in Fig.l, where the regime which we analyse is the upper half- 
plane (U>0). The superconducting (SC) instability which is present in this diagram 
is ignored in our RPA approach. However, since we are interested in the the effects of 
magnetic field, this omission can be justified even in the case when for = the singlet 
SC state (SS in Fig. 1) is stable, i.e. when z/ < —1/3. Namely, the superconducting 
phase is suppressed by the orbital effect of a magnetic field. The critical field at which 
the critical temperature for the singlet SC state drops to zero is given by ||2^ 

Hf = , (15) 

7\/2C(3)/iBi;it 

where Tsc is the critical temperature for the singlet superconducting state in the absence 
of magnetic field. Considering {gi + g2) as the corresponding effective coupling constant 
|jl|], one easily finds that Tsc is related to the critical temperature for the charge density 
wave at zero magnetic field, T°, by 

Ttvp 1 — 3z/" 



Tsc = T° exp 



U l + 3z/ 



(16) 



Equations (|1^, |16D give the estimation for the magnetic field above which our RPA 
results are valid even in the regime when the singlet superconductivity overwhelms the 
CDW. 

It is useful for further discussion to mention here a few characteristic possibilities 
regarding the value of the parameter v in the CDW (i.e. [/ > 0) systems. Taking into 
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acount only pure backward electron-phonon interaction one has z/ = 0. The inclusion 
of the presumably weaker repulsive Coulomb interaction between electrons shifts v to 
some positive value. From the other side, a pure Hubbard model with attractive on-site 
interaction corresponds to v = —1. Altogether, u covers a wide range of theoretically 
allowed values, but it should be noted that in the most frequent electron-phonon CDW 
systems this range is limited to z/ > 0. 

Finally, it should be noted that the function /H from the matrix (^) can be expressed 
in the factorized form, /H = Thus, for f/ > the condition ( ]TT| ) reduces to 

X:'(q,Te) = 0, (17) 

i.e. to the divergence of the susceptibility (M_m1). Indeed, in the limit H — >^ the 
component M_ reduces to the pure CDW component M4 and the divergence of X- 
coincides with the condition for the CDW instability, 1 + f/cXo = 0. Since the further 
discussion involves only the ordering with finite components M3 and M4, we simplify 
the notation for TJ' and qy in eq.(]T7|) by skipping the index ||. 



3 Discussion 



For = 0, the wave vector oi CDW ordering for H = is defined by the maximum 
of the susceptibility XoiSii = 0) (H) in the limit qo 0. Of course it is located at q = 
0, i. e. at the wave vector of perfect nesting. The corresponding critical temperature 
T° = (27i?i?/7r) exp{—nvF/U) defines the temperature scale of the problem. 
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We want now to calculate the position of the minimum of xl^ (eq.|T3|) in the mo- 
mentum space. The criterion of local stability of the ordering with q = at finite H 
can be derived from the quadratic expansion of with respect to and qy. For 
T = Tc{H) it suffices to expand the [..] bracket in eq. (p!3|) . Noting that there is no 
bilinearly mixed term (qxqy), one gets 

U~\gXZ' ^ -Xo + a^ql + ayq^ + + 0{qlql q^, ..) , (18) 



with 



2xr {i-^wxi 

2 dq: 



'dxoV 
, dqx I 



(19) 



= -ol^' (20) 



and 



1^0 
4! dql 



h = -771^ ' (21) 



with the values of Xo = Xoi^lx, %, = 0) and its derivatives taken at q^ = 2qp, qy = 0. 
For the later purposes we include one (~ q^) of fourth order terms in the expansion 



(|1|). Note that the expansion (|T8D is valid for q^ < AnT^/vp^qy < AnT^ / (Ub) . 

The dependence of the critical temperature for the ordering at q = on the mag- 
netic field follows from the equation 

Uxo = 1. (22) 

For small values of H this expression reduces to the known result for the suppression 
of the critical temperature due to the Pauli splitting of the electron band 

T, = TO (l - 7C(3)(/isi7/27rT°)2) . (23) 
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The dependence of the coefficients ax, ay, and by on the magnetic field follows 
straightforwardly from eqs.(||, |^). To this end we use the relation 



-KVp 



(24) 



where ^ denotes the digamma function, and expand the coefficients (|^) (with = 0) 



in terms of qy up to the quartic contribution. Taking into account also eq. (|22|) one gets 
at T = T^{H) 

1 



with VI/ ^.1/(1 + ^), 



and 



•9' o T 2z/[/ ( d ^ \- 
dq^ (1 - u)txvf \oqx J 



qx=2qp 



(25) 



1 ( h V,2 



(26) 



^ Tivp \2nTcr]h J 
The coefficients ay and Py in eq . (|27|) are given by 



(27) 



ay = Re"^ (l + ih{l+7]/2)) + Re^ (l + ih{l-7]/2) 



-2Re^! i^ + ih 



{21 



and 



Py ^ li?eM/Q + ^/i(l + r/)) + ii?ev^Q + z/.(l-r^) 

-^Re^H Q + ^/^(l + r//2)) - i/^ev]/ Q + zh{l - r,/2) 



3 /I 

+2** 2 



(29) 
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respectively. Here we have introduced the dimensionless variable h = ^bH /2'kT . Note 
that the quantities and ay determine the longitudinal and transverse correlation 
lengths for CDW fluctuations = U^/a^ and = U ^Ja^ when the temperature is 
close to Tf.. 

The wave vector of ordering stays at q = as far as the coefficients a^iji) and 
are positive, and starts to move in the longitudinal or transverse direction when the 
former or later coefficient changes sign. As it is seen from eqs. (|25|p6D , the function 
ax{h)/ax{h = 0) contains the interaction parameters U/vp and z/, while the function 
ay{h)/ay{h = 0) depends only on the ratio r] = qo/qp = ebvp cos9/ fiB which mea- 
sures the relative impact of the orbital and Pauli coupling on the CDW. Note that 
the parameter rj can be easily changed by varying the angle 6 between the direction 
of magnetic field and the c axis. Since the reasons for possible deviations of stable 
components q^ and qy from zero are essentially different, it is appropriate to consider 
each case separately. 

The longitudinal component of the CDW wave vector. The coefficient a^ changes 
its sign at the critical field h^x = ^sHcx/'^'^T shown in Fig.2. For small values of v 
this dependence is given by 



with /i° = ^bH°/ {2t{T) = 0.304. As it is seen in Fig.2, all curves hcx{v) pass through 
two common points, given by z/ = 0, h^x = 0.304 and u = 1, hex = 0. At these points, 
hex does not depend on U. The first point (z/ = 0) corresponds to the CDW ordering 




(30) 
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with the SDW couphng equal to zero. In the second case (z/ = 1) the interactions in 
the CDW and SDW channels are of equal strengths and opposite signs, i.e. we are 
at the CDW — SDW boundary. There, the longitudinal splitting of the wave vector 
starts already at h^x = 0. The change of sign of at hcxi^i U) causes a second order 
transition from the phase with q=0 named CDWo, to a phase with q=(g^(/i),0), named 
CDWx- The dependence of the wave vector on the magnetic field is the solution of the 
equation dxl^ /dqx = 0, and can be written in the form 

Qx = ^^fu,uih). (31) 

Vf 

The function fu,u{h) is shown in Fig.3 for U/ttvf = 0.2 and few values of u. In the 
limit h » 1 one has qx{h) 2^bH/vf, as it was already shown previously in the 
case of repulsive Hubbard model (z/ = — 1) [^]. 

The transverse component of the CDW wave vector. The critical field hey at which 
a finite transverse component of the CDW vector develops is shown in Fig.4. The line 
hcy{f]) corresponds to the second order transition from CDWo to a phase CDWy with 
a transversely shifted wave vector. At small values of the dependence hcy{ri) is given 
by 



hey ^ K^l + 0.088r/2. (32) 

Here we use the approximative expression Re'^(^ + ix^ ~ ^ (i) ~'~ 8-414x2(1 + 
3.81a;^)~^, valid for small values of the argument x. The dependence of the wave 
vector component Qy on h for a fixed value of rj can be represented by 



2 

Qy = - arcsm 



tb 



(33) 
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with the function (7^ (h) shown in Fig. 5. Note that unhke hey the wave vector component 



Qy depends on ti,. For small values of h — hcy the function g^^ reduces, after using eq.(^ 
to 



g,{h)^2r^h^-^^, (34) 

with ay and Py given by eqs.(|^) and (|29|) respectively. On the other hand in the 
high field limit h ^ h^y and for rj = the function gri=o{h) is asymptotically given by 
gr)=o{h ^ 1) — s> /i + K, where k is of the order l/vr. Note that the transverse shift of 
the wave vector does not depend on the interaction (z/ or U). It however depends on 
1], i. e. on the relative impact of the Pauli and orbital effects. The reason for this is in 
the fact that all interaction dependence enter with 6 (see eq.|^), which is equal to zero 
if = and if the nesting is perfect. Thus, only a phase CDW^ is affected by the 
finiteness of the SDW coupling constant Us (i. e. z/). 

The phase diagram for h larger than hex and/or hey To provide some ideas on 
the variation of the wave vector of instability at magnetic fields larger than critical 
values hex and hey, it is useful to consider the symmetry and the shape of functions 
Xo(q),Xg(q) and xl^(q) at strong magnetic fields {h of order 1). Note firstly that all 
these functions are even in Qx and Qy. The function Xo ^(^) at T << tf, has the line of 
local maxima given by pl[ 



qx = ± 



4h . %b I (T 
— sm — — — O — 

V 2 V \2tb 



(35) 



When the Pauli term is introduced, the maxima of Xt(^) ^ill move to the left and 
those of Xj(q) to the right by 2qp along the axis g^. The lines of local maxima of the 
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susceptibilities Xo, X? Xi ^i'g shown in Fig.6a. For h large enough the function 
Xg(ct) = ^/X^Xl: together with the function — xl^(q), will have two pairs of degenerate 
maxima in q-space as candidates for absolute maxima. These two pairs have approx- 
imate positions at (±2gp,0) and (0, ±| arcsin(fpgp/2tfc)) (denoted as A,A' and B,B' 



respectively in Fig. 6a), in accord with the asymptotic limits given by eqs. (pT] , |33D . In 
Fig. 6, we also show the function — xl^(q) for three characteristic choices of parame- 
ters u, 1] and U, i.e. when the absolute maxima are at {0,±qy) (Fig.6b), at (±2gp,0) 
(Fig. 6c), and when the two pairs of maxima have the same value (Fig.6d). As it was 
shown above, the phase transitions from the CDWo to CDWx and CDWy (Figs. (6b) 
and (6c)) are of the second order. The transition between the orderings CDW^ and 
CDWy, caused by the competition of two maxima in — xl^(q) (Fig.6d) is of the first 
order since the wave vector has a discontinuous jump between points (g^.', 0) and (0, qy) 
(i.e. between points A and B in Fig.6). 

To complete the phase diagram it is necessary to calculate the magnetic field de- 
pendence of the critical temperatures, defined as the solutions of eq.(|) for q = (0, 0), 
q = (g^, 0) and q = (, qy), and denoted by Tq, and Ty respectively, and to determine 
max\TQ{H),Tx{H),Ty{H)]. The dependence of critical temperatures Tq^T^ and Ty on 
H for few values of v and r] and for Uj-nvp = 0.2 is shown in Fig. 7a. The sections 
of lines Tx{H) and Ty{H) determine the critical magnetic fields and the temperatures 
of the first order transitions. Note that the present analysis is based on the Landau 
expansion 

F = Jd'q x:^(q)[M_(q)M!(q)] + 0{{Mt}), (36) 



which is restricted to the range of temperatures not far below 



Since the complete phase diagram depends on three parameters, H, v and 77 (with 
fixed f/), it is appropriate to use two planes, (if, z/) (77 being a parameter) and (if, ri) [y 
being a parameter), for its presentation, as shown in Figs. (7b) and (7c) respectively. 
We stress a particularly interesting situation v ^ \~ for which the critical field hex 
goes to zero, and three phases [CDWq^ CDW^ and SDW) are present in the narrow 
range of parameter u. Note also the presence of the point in Figs. (7b) and (7c) at which 
the CDWo, CDWx and CDWy orders meet. The dependence ulrj) which defines this 
tricritical point is shown in Fig.(7d). The corresponding magnetic field weakly varies 
with u (i. e. rj), as is seen in Figs. (7b) and (7c). The line in Fig.(7d) thus divides 
the region where the wave vector shifts firstly in transversal direction from the region 
in which only a longitudinal shift is possible. Furthermore, among the CDW phases 
from Figs. (7b, c) only the phase CDWx has a finite fraction of the component M3 [see 
eq.(|T2|)], and is thus a CDW — SDW hybrid. The ratio of components M3 and M4 
follows from the constraint M_|_ = 0. At T = T^x it is given by 



and shown in Fig. 8 for few values of z/. Note that IM3/M4I tends to 1 as one approaches 
the CDWx - SDW transition. 

Influence of the imperfect nesting. Let us finally consider a case when the imperfect 
nesting is introduced through a finite effective next-nearest neighbor hopping t'j,, which 



max[T^{H),T^{H),Ty{H)]. 




(37) 
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can be usually increased by e. g. applying a strong pressure on a CDW system. For 



example the relevant pressure scale in NhSe^ is about 10 kBar|jT2|, ^ . 

At small values of the magnetic field, the the critical temperature for the phase 
CDWq can be readily found from the eg. (^2]) yielding 

'dxo 



T — T° ~ — 



dT 



2 dql 



(3^ 



The values of ^ and as functions of t'^ are given in ref. [0. For small t'j,, the 
coefficient is given by 



vf 

a.. ' 



327r3T° 



-^"(1/2) + ^'^(1/2) 



27rT° 



(39) 



where ~ —16.83 and —771 AT. As one sees from the eg. (^) , the orbital 

and Pauli effects are in competition, the former trying to enhance, and the latter to 
suppress T^. For small t'^/T°{t'f, = 0) the function is proportional to t'^^. Moreover, 
the imperfect nesting decreases the coefficient a^. Altogether, the general trend of the 
small t'f^ is to flatten the H dependence of Tc. 

For the sake of space, we present the result for the critical temperature wich follows 
from the eg. ([Tl|) only for the case of attractive Hubbard interaction (i.e. U/tcvf = 0.2; 
u = — 1)|2^ and for t] = 2.5. In this regime the orbital effects are strong enough, 
which excludes the stabilization of the CDWy ordering when the nesting is good. The 
interplay between two effects of a magnetic fleld is a main characteristic of the phase 
diagram for imperfect nesting, given in Fig.9. As the parameter t'j, increases from 
zero the critical temperature only monotonously shifts to lower temperatures, still 
decreasing with a magnetic field. In other words, our results for the perfect nesting 
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can be applied even to the systems with a moderate finite imperfect nesting, i.e. when 
the critical temperature remains far above the value of t^. The orbital effects enter 
manifestly into play at rather large values of t'^, enhancing the critical temperature 



initially, as it was observed in NbSe-^ |T^. The eventual supression of by Pauli 



term at high magnetic fields will make these diagrams basically different from the 
mean- field one for the FISDW with the orbital coupling only |17[, where no eventual 
supression of the is present. For a very bad nesting, i. e. for t'^ comparable to t'^ 
[where t'^ ~ T°{t'f^ = 0)) is the imperferct nesting parameter at which the CDW is 



destroyed at zero field fT^]]; the cyclotron frequency becomes the first relevant energy 
scale, giving the rise to a cascade-like shape, associated with the quantized field induced 
CDW phases. Notice that our approach does not explain the strong field breakdown 
of the high field phase in {TMTSF)2C104 p2|, since the Pauli term does not affect the 
SDW. 



4 Conclusion 

The main result of the present work concerns the phase diagram of a CDW system 
in an external magnetic field. There are three physical parameters which characterize 
this diagram, namely the ratio of the SDW and CDW coupling constants, the strength 
of the magnetic field and its direction with respect to the most conducting plane (x, y). 
The respective parameters are u, h and t]. We recall that rj also measures the relative 
impact of the orbital coupling with respect to the Pauli coupling. 
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In the case of good nested Fermi surface the wave vector of the CDW has a general 
tendency to shift from its zero field value {2kp,7i/b) as the magnetic field increases 
[see Figs.(7a-c)]. This shift starts continuously, and may occur either longitudinally 
or transversally with respect to the chain direction. The longitudinal shift is governed 
solely by the Pauli coupling, with the corresponding CDWx state being a hybrid of 
the pure CDW and of the SDW component parallel to the magnetic field. Both, the 
critical value of the magnetic field hex at which starts to shift, and the relative weights 
of the CDW and the SDW, depend on the ratio u. Both, and the CDW-SDW 

hybridization increase with the magnetic field. It is important to mention that hex, 
Qx and the hybridization ratio do not depend on U because all mean field properties 
concerning a longitudinal tilt of the wave vector are given by pure one- dimensional 
expressions. 

The shift of the CDW wave vector in the transverse direction is affected by both 
orbital and Pauh couplings. Contrary to the CDWx, the CDWy is not a CDW-SDW 
hybrid, and therefore is not influenced by the parameter v. It exists only when is 
finite, although the critical magnetic field h^y does not depend on t;,. However, tb 
influences the variation of qy aX h > hey, as shown by eg. (P3|) . qy{h) decreases with tb 
and increases with the magnetic field. According to the general fact that the orbital 
effects lower the dimensionality of the electronic motion ||16|, the effect of the increasing 



T] is to favor the CDWx- After some critical value of rj (dependent on z/), the orbital 
impact reduces the phase diagram to the pure one- dimensional one, consisting only of 
the CDWo and CDWx, as it is seen from Fig. 7b. 
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At = and for z/ < 0, the shift of the wave vector is at first directed perpendicu- 
larly, and jumps to the longitudinal direction at some higher magnetic field. This jump 
between CDWy and CDW^ is a first order transition. On the contrary, for < z/ < 1, 
the wave vector is shifted longitudinally for all magnetic fields higher than the criti- 
cal field hcx{t^)- Furthermore, hcx{i^) tends to zero as u approaches unity. The point 
if = 0, z/ = 1 is therefore tricritical, since z/ > 1 is the range of SDW stability. 

The Pauli and orbital terms together cause a rather complex magnetic field depen- 
dence of the critical temperature in systems with a finite imperfect nesting. This is 
illustrated in Fig. 9 in which t'f^ varies from zero to the range above the critical value 
4*, at which the CDW ordering is completely eliminated at zero magnetic field. A 
rich dependence Tc{H) contains the suppression by the Pauli term, enhancement by 
the orbital effects and, for large values of t'^, a cascade-like shape characterizing the 
field-induced DW. This phase diagram is quite general and not limited to the value 
V = —1, chosen in Fig. (9). 

Our analysis for the perfect nesting case, showing a strong dependence of the critical 
properties in magnetic field with the ratio i>, could find an appropriate experimental 
support e.g. in the MX compounds. The low- dimensional nature of these materials 
corresponds to our model. From our analysis, a particularly interesting possibility is 
that the Coulomb and electron-phonon forces can be tuned in a predictable manner 
by external pressure or chemically |2^, allowing to approach the phase boundary 



between CDW and SDW, corresponding to u = 1. As we approach the boundary from 
the CDW side, the critical field for the CDWo CDW^ transition hex = Hcx/2nT 
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will decrease rapidly toward zero, regardless of the value of T ^ Tc- Even for large Tc, 
by adjusting carefully Hex can decrease to experimentally reachable values, being 
extremely sensible to the variation of the parameter v. We point out that a search for a 
magnetic field induced phase transitions in a CDW phase with strong SDW fluctuations 
(introduced by high pressure, for example) could confirm our predictions. 

In NbSe^ a phase transition in the 59K CDW phase induced by magnetic field was 
found by observing that a threshold electric field for the collective CDW motion 
is strongly reduced when magnetic field increases beyond the critical point. The naive 
explanation that this is a simple CDWo CDWx one-dimension-like transition due 
to only Pauli term must be taken with caution. Namely, the observed effect strongly 
depends on the angle 9, indicating that the orbital effects are also involved. This might 
mean that the strong orbital contributions, provided by a badly nested Fermi surface, 
affect the phase diagram. However, we believe that the Pauli term has an important role 
in this transition, since it enables the shift of the wave vector from its commensurable, 
perfect nesting position. We remind that pure orbital effects can affect the wave vector 
only if it is not at the perfect nesting position [like in e.g. {T MT S F)2C10 [P^]. 
The fact that the nesting in NbSe^ is quite bad can be deduced from relatively strong 
pressure dependence of Tc {dTc/dP ^ —6.25K/Kbr)^^. Indeed, from the comparisson 
of a very weak enhancement of Tc with magnetic field with our results in Fig. 9, it 
follows that the value of should be rather large. 

Finally, our analysis of the imperfect nesting case can somewhat enlighten the re- 
cent measurements in the compound Per2[Au{mnt)2], where the suppression of the 
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critical temperature proportional to the square of the magnetic field was found, but 
with a coefficient smaller than that which follows after taking only the Pauli coupling 
and a perfectly nested Fermi surface f^. From eq.(^) and from Fig. 9 one can conclude 
that the reason for the fiattening of the suppression of is just the finiteness of t'^,. 
However, the situation is not so simple. At finite values of t'^, the orbital effects come 
into play, in contrast to the experimental results which are independent on the field 
direction. If we just ignore the orbital effects, we get t'^ ~ 7AK as a imperfect nesting 
parameter fitting the experimental curve. Finally, we indicate that the measurements 
of the critical properties in a magnetic field, and with pressure large enough to al- 
most or completely destroy the zero-field CDW, could show very strong, cascade-like 
enhancement of the Tc for the quantized field induced CDW phases. 
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FIGURES 

FIG. 1. A phase diagram of the one-dimensional system in {v-iU) plane, in the 
absence of a magnetic field. 

FIG. 2. Scaled critical magnetic field hex = f^sHcx/i'^'^T) as a function of the 
parameter u for few choices of the coupling constant: U/ttvf =0; 0.2; 0.4; 0.6. Note 
that one has to insert T^^H) [and not Tc{H = 0)] into the defining expression for hex 
in order to get a phase diagram with H dependence. 

FIG. 3. The function fu,u{h), determining the dependence of the longitudinal shift 



of the wave vector on the magnetic field (see eq.31), for U/nvp = 0.2 and v = —1 (A), 
-0.5 (B), 0. (C), 0.25 (D), 0.5 (E), 0.75 (F), 0.99 (G). 

FIG. 4. Scaled critical magnetic field hey = f^BHcy/ {2nT) as a function of the 
parameter rj. 

FIG. 5. A function gr]{h), determining the dependence of the transverse shift of the 
wave vector on the magnetic field (see eqJS^), for r/ = (A), 0.5 (B), 1 (C), 1.5 (D), 2 
(E). The inset shows the large x behavior for a case with no orbital effects [rj = 0). 

FIG. 6. a) The lines of local maxima in q-space of the susceptibility Xoilx, %) (full 
line) without magnetic field, of xt = Xo{<ix + 2gp, qy) (dot-dashed line) and of Xi = 
Xoilx—'^lp, %) (dashed line). (A,A') and (B,B') are the two pairs of degenerate maxima 
of Xg{.1x,(ly)- The figures b), c) and d) show the function — xl^(q) at T = 0.42T° and 
^ibH = 1.14T° respectively for three cases: u = —1 {CDWy is stable, provided that 
the maxima are in the points B,B'), u = —0.1 {CDWx is stable; the points A,A' are 
dominant) and u = —0.33 (the first order critical point between CDWy and CDWx, 
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since all four points, A,A',B and B' are of the equal height). 

FIG. 7. a) The critical temperatures Tq, and Ty for the CDW instabilities with 
the wave vectors q= (full line), q=(?a;,0) (dotted line) and q = (O^Qy) (dot-dashed 
line) respectively. For all curves U/irvp = 0.2. 

b) The phase diagram in [v, jiBH / {2t:T°)) plane for r] — The changes of the 
diagram with finite r] (here, r] — 1) are shown by dashed lines. 

c) The phase diagram in {rj, hBH/{2ttT^)) plane for u = —1. The changes of the 
diagram when u — —0.5 are dashed. 

d) The curve in the {i',')]) plane which defines the tricritical point in Figs.(b) and 

(c). 

FIG. 8. The relative weight of the SDW and CDW components in the hybrid 

phase CDWx as a function of magnetic field, for U/tivp = 0.2 and for u = —1 (A), 
(B), 0.99 (C). 

FIG. 9. The critical temperature vs magnetic field for a series of values of tyt'^ 
and for U/nvp — 0.2 and u — —1. 
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